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Abstract. We give a more geometric approach to an algorithm 
for deciding whether two hyperbolic 3-manifolds are homcomor- 
phic. Wc also give a more algebraic approach to the homeomor- 
phism problem for geometric, but non-hyperbolic, 3-manifolds. 

1st November 

1. Introduction 

The homeomorphism problem for closed orientable triangulated 3- 
manifolds has been studied for many years, with partial results by 
many authors. The work of Perelman [bll [T5| [16] proving Thurston's 
Geometrization Conjecture finally allowed a complete solution for irre- 
ducible such manifolds, which has been described by Jaco in his Beijing 
Lectures [7] , and also in the first chapter of the book by Bessieres et al 
[3]. The work of many previous authors is put together, and different 
algorithms are used to deal with the Haken case, to find the JSJ de- 
composition, and to deal with different geometries. In the case of two 
closed hyperbolic 3-manifolds Mi and M 2 , Mostow rigidity tells us that 
Mi and M 2 are homeomorphic if and only if their fundamental groups 
are isomorphic. Thus the homeomorphism problem for Mi and M 2 
can be solved by appealing to Sela's solution [19J of the isomorphism 
problem for torsion-free word hyperbolic groups. The initial aim of this 
paper was to give a more geometric approach to the homeomorphism 
problem in this case, which avoids quoting Sela's work. But in addi- 
tion, we also give a more algebraic approach to some other parts of the 
homeomorphism problem. We are not claiming that the algorithms 
we present are in any way superior to those referred to by Jaco and 
by Bessieres et al. Our aim is simply to increase the range of appli- 
cable algorithms from which to choose. In this paper, we will mostly 
consider closed, orientable and irreducible 3-manifolds. It is known 
that given two triangulated closed orientable 3-manifolds Mi and M 2 , 
there is an algorithm to find the geometric structures on the geometric 
pieces of each, and then there is an algorithm to decide whether or 
not the pieces of Mi are homeomorphic to the pieces of M 2 (see e.g. 
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[TTj 1.4]). This algorithm is described in [3J 1.4.1], where it is claimed 
that this solves the homeomorphism problem for triangulated closed 
orientable 3-manifolds. (A similar claim is made by Matveev in his 
book [T2"| p. 214], using his method of spines for Haken manifolds.) It 
was pointed out to us by Henry Wilton, that there remains an orien- 
tation problem when considering connected sums. If M and N denote 
two closed orientable 3-manifolds, and N denotes N with the opposite 
orientation, it is not clear how to decide whether the connected sum 
M#N is homeomorphic to M#N (though the geometric pieces are 
clearly homeomorphic). 

In [TT] , Manning gave an algorithm to decide if there exists a hyper- 
bolic structure on a closed orientable 3-manifold given by a triangula- 
tion. If there is such a structure, then Manning's algorithm constructs 
a finite sided polyhedral fundamental region in hyperbolic 3-space. In 
the first section of this paper we use Manning's work to give a new 
algorithm to decide whether or not two closed 3-manifolds Mi and M<i 
are homeomorphic, when they are given by triangulations and known 
to be hyperbolic. Manning's algorithm will construct a finite sided 
polyhedral fundamental region in H 3 for each manifold Mj, and we 
show how to estimate how many copies of P2 must be glued together 
to contain Pi and vice versa. This allows one to bound the number of 
maps from the generators of tti(Mi) to 7Ti(M2) which might give rise to 
an isomorphism, and vice versa, enabling one to check each such map 
in turn and decide whether or not it is an isomorphism. 

In the second section of this paper, we take a more algebraic point of 
view. The main ingredients are the existence of a bi-automatic struc- 
ture on the fundamental groups of most geometric 3-manifolds and 
the fact that Perelman's work essentially reduces the homeomorphism 
problem for closed orientable irreducible 3-manifolds to the isomor- 
phism problem for their fundamental groups (with the notable excep- 
tion of lens spaces). 

Given a triangulated closed orientable 3-manifold, Jaco and Tollef- 
son [9j §7] give an algorithm producing a decomposition into irreducible 
pieces. When the manifold is Haken, they give [9j §8] algorithms for 
finding its JSJ decomposition and finding the Seifert invariants for 
each Seifert piece. Thus for a Haken closed Seifert fibre space, their 
algorithm will yield all the Seifert invariants and hence a geometric 
structure. Now suppose M is a closed orientable 3-manifold which is 
known to satisfy Thurston's Geometrization Conjecture. In [TTJ 1.4], 
Manning combined the algorithms in [9] with his algorithm for finding 
hyperbolic structures to obtain an algorithm which finds the geometric 
structures on all the pieces of M. In this section we describe a some- 
what simpler algorithm, with substantial algebraic ingredients, which 
can answer a somewhat simpler question. Namely, if M is a closed 
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orientable 3-manifold which is known to be geometric, then our algo- 
rithm can decide on which geometry M is modelled. To answer this 
simpler question, we do not need Manning's algorithm for hyperbolic 
manifolds, and our procedure differs from that of Manning in the other 
cases. One natural application of our algorithm would be when M is 
orientable and irreducible and not Haken, as then Perelman's work im- 
plies that M must be geometric. Our approach also works in the Haken 
case and is rather different from that in [9]. Our algorithm starts by 
finding a bi-automatic structure for the fundamental group of M, when 
such a structure exists. As for many algorithms for 3-manifolds, the 
main topological ingredient is an algorithm due to Jaco and Oertel [5] 
to decide whether a triangulated 3-manifold is Haken, and if so, to 
produce embedded incompressible surfaces. 

2. A NEW ALGORITHM TO DECIDE IF TWO HYPERBOLIC 
3-MANIFOLDS ARE HOMEOMORPHIC 

Let Mi and M 2 be closed 3-manifolds each given by a triangulation 
and known to be hyperbolic. It is straightforward to use these triangu- 
lations to write down finite presentations for 7r 1 (M 1 ) and 7r 1 (M 2 ). In 
order to apply Manning's algorithm in [TTJ, we first need an algorithm 
for solving the word problem in each 7i"i(Mj). As Mj is known to be hy- 
perbolic, its fundamental group is word hyperbolic and so bi-automatic. 
Then, the algorithm of jH 3.4.1] (extended to the bi-automatic case, 
as stated in [5] after Lemma 8.2) produces a bi-automatic structure. 
Alternatively Papasoglu [13] gives an algorithm for calculating the hy- 
perbolicity constant 5 of the presentation, and with that a solution to 
the word problem is easily built as in [TJ. It follows that there is an 
algorithm which can be applied to the given presentation for 7Ti(Mj) 
which will find the bi-automatic structure. (We note that chapter 5 
of [I] gives the algorithm to find automatic structures; the additional 
languages and axioms needed to extend to bi-automatic structures are 
easily added.) Once this is done, there is an algorithm for solving the 
word problem in 7Ti(Mj) (alternatively one can use the fact that the 
groups are residually finite). Now, for each Mi, Manning's algorithm 
in [11] constructs a convex finite sided polyhedral fundamental region 
Pj in H 3 . His algorithm describes the vertices, edges and faces of Pj. In 
addition, it describes the face pairing isometries needed to recover Mj 
from Pj. It is now straightforward to write down a new presentation 
for TTi(Mi) with the face pairing isometries of Pj as generators. In what 
follows we will use these presentations rather than the ones obtained 
from the initial triangulations. 

Now we consider the tiling of H 3 by translates of P%. Given a union 
X of such translates, we let star(X) denote the union of all translates 
of P2 which meet X in at least one point, and let star n (X) denote the 
result of applying the star operation n times to X. Thus star 2 (X) = 
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star (star (X)). It may sometimes be convenient to write star°(X) = 
X. We consider the sequence star n (P 2 ), n > 0, of subsets of H 3 . As the 
union of these subsets equals H 3 , their diameters must tend to infinity. 
We will refine this obvious fact in the following way. We will show in 
Lemma 12.11 below that there is an algorithm to find a positive number 
R, such that star n (P 2 ) contains the metric ball B(P2,nR). 

Assuming this lemma for the moment, we now proceed as follows. 
If Mi and M 2 are homeomorphic, Mostow's Rigidity Theorem implies 
that they must be isometric with their hyperbolic metrics. This yields 
an isometry ip from H 3 tiled by translates of Pi to H 3 tiled by translates 
of P 2 . By composing with the action of an element of 7Ti(M 2 ) on H 3 , 
we can suppose that ip(Pi) meets P 2 . Let d\ denote the diameter of 
P\. Then tp(Pi) must be contained in the metric ball B(P 2 ,di). If n 
is an integer such that nR > di, it follows that (p{Pi) is contained in 
star n (P 2 ). Let a be a face pairing isometry of P\. Then the isometry 
(paip^ 1 of H 3 pairs faces of <p{Pi). As a lies in 7Ti(Mi), the isometry 
ipoap~ l lies in tti(M 2 ), and so preserves the tiling by translates of P 2 . 
Hence (poop -1 must send a certain translate of P 2 which is contained 
in star n (P 2 ) to another such. Pick a path in star n {P 2 ) which joins the 
interiors of these two translates of P 2 , does not meet any edges and is 
transverse to the faces, and let N denote the number of times this path 
meets a face. Then (poop -1 can be written as a word in the face pairing 
generators of tti(M 2 ) of length N. Let k\ denote the maximum number 
of translates of P 2 around an edge, and k 2 denote the maximum number 
of translates of P 2 around a vertex, and let k denote the maximum of 
k\ and k 2 . Then in star(P 2 ), we can join P 2 to any translate of P 2 , 
by a path which crosses at most k faces. Note that a translate of P 2 
in star(P 2 ) may have just a vertex in common with P 2 . It follows 
immediately that in star n (P 2 ), we can join P 2 to any translate of P 2 , 
by a path which crosses at most nk faces. Hence, in star n (P 2 ), we can 
join any two translates of P 2 , by a path which crosses at most Ink 
faces. We conclude that if there is an isomorphism from 7Ti(Mi) to 
tti(M 2 ), there is one which maps each face pairing generator of 7Ti(Mi) 
to a word of length no more than Ink in the face pairing generators of 
7Ti(M 2 ). Similarly we can find integers m and k' such that if there is 
an isomorphism from iri(M 2 ) to 7Ti(Mi), there is one which is inverse 
to the previous one, which maps each face pairing generator of ir 1 (M 2 ) 
to a word of length no more than 2mk' in the face pairing generators 
of 7Ti(Mi). This gives us a finite list of possible maps from generators 
of 7Ti(Mi) to elements of tti(M 2 ), and vice versa. For each such map 
we can check whether it is a homomorphism, and for each pair of such 
maps, can check if their composite is the identity. Each of these checks 
again requires the solution of the word problem. Thus we can check 
whether or not there is an isomorphism between 7Ti(Mi) and iri(M 2 ). 
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Before giving the proof of Lemma 12.14 we need to discuss algorithms 
for bounding various distances related to the convex polyhedra Pj. 

We will consider the upper half space model of H 3 , and all coordi- 
nates used will be euclidean. Recall that in this model a hyperbolic 
geodesic is either a vertical line, or a semi-circle in a vertical plane 
centred at some point of the base plane M 2 of the model. Also a hy- 
perbolic plane in this model is either a vertical plane or is a euclidean 
hemisphere centred at some point of the base plane M 2 of the model. 
The two points which form the intersection of such a semi-circle with 
the base plane M 2 of the model will be called the boundary of the geo- 
desic, and the circle which is the intersection of such a hemisphere with 
the base plane M 2 of the model will be called the boundary of the hy- 
perbolic plane. If a hyperbolic geodesic is a vertical line in this model, 
its boundary consists of one point in M 2 and one point at infinity. If a 
hyperbolic plane is a vertical plane in this model, its boundary consists 
of a line in M 2 and one point at infinity. 

We recall that Manning's paper [11] produces the hemispheres which 
contain the faces of Pi, and that each of these hemispheres has centre 
with coordinates which are algebraic numbers and has euclidean radius 
which is also an algebraic number. If the intersection of two of these 
hemispheres is non-empty, it is a semi-circle whose boundary points 
have coordinates which are algebraic numbers and whose euclidean ra- 
dius is also an algebraic number. Further the vertices of Pi, each of 
which is the intersection of three of these hemispheres, also have coor- 
dinates which are algebraic numbers. These numbers can be approx- 
imated to any required degree of accuracy over the rational numbers 
(to lie within an "isolating interval" with rational endpoints) using 
standard methods of symbolic computation, as described by Manning 
(with reference to [2] and [ID])- The algorithmic process starts from 
the results of Manning's algorithm, which gives coordinates (i.e. ratio- 
nal estimates and isolating intervals) for the vertices of the polyhedra 
Pi, and for the centres of the semi-circles and hemispheres defining the 
edges and faces of Pi, and estimates of their radii. Suppose that all 
these isolating intervals are of width at most e, and to begin, suppose 
that e has been chosen so that 1/2 C+1 < e < 1/2 C . We refer to this 
number as the error in our calculations. 

Recall that we work in the upper halfspace model for HI 3 , and we 
will use euclidean geometry and euclidean coordinates, radii etc. After 
estimating euclidean distances, we can then estimate hyperbolic dis- 
tances using the usual formulae. For simplicity in the following seven 
statements, we will say that a point in the upper half space model of 
HI 3 is algebraic if its coordinates are algebraic, that an infinite geodesic 
in this model is algebraic if its boundary points are algebraic (in the 
base plane z — 0), that a hemisphere in this model is algebraic if its 
euclidean centre (in the base plane z — 0) and radius are algebraic, and 
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that a vertical plane is algebraic if it contains three algebraic points 
(or has two algebraic boundary points). Finally a compact geodesic 
segment contained in an algebraic geodesic is itself algebraic if its end- 
points are algebraic. 



(1) The distance between two distinct algebraic points is not algo- 
rithmically computable, but we can compute a positive lower 
bound which will suffice for our requirements. The error in- 
curred here is at most 2e. 

(2) Finding the distance between an algebraic point X and a dis- 
joint algebraic geodesic A in the hyperbolic plane can be algo- 
rithmically reduced to finding the distance between two alge- 
braic points as follows. Let \x denote the semi-circle through X 
which meets A at right angles. We can find the centre of this 
semi-circle on the base line R by solving quadratic equations 
specifying that the centre is equidistant from X and the point 
A fl yU, and that the line from the centre to A D \i is tangent to A. 
Hence we can also find Afl/x. Now the required distance equals 
the distance from X to the point A n fi, so we can apply 1). 

(3) In the same way, finding the distance between two disjoint al- 
gebraic geodesies (without a common boundary point) in the 
hyperbolic plane can be algorithmically reduced to 1). 

(4) In the same way, finding the distance between an algebraic point 
X and a disjoint algebraic hemisphere (or vertical plane) II in 
hyperbolic 3-space can be algorithmically reduced to 1). 

(5) In the same way, finding the distance between two disjoint al- 
gebraic hemispheres (or vertical planes) whose boundaries are 
also disjoint can be algorithmically reduced to 1). 

(6) To find a lower bound for the distance between an algebraic 
semi-circle A and a disjoint algebraic hemisphere II, such that 
the boundaries of A and IT are also disjoint, we will choose an 
algebraic hemisphere E which contains A, and is disjoint from 
II, so that their boundaries are also disjoint. Once such E has 
been found, the distance between E and IT, which can be found 
as in 5) gives the required lower bound. 

We can choose E as follows: let P and Q denote the boundary 
points of A in the plane z — 0, and let L be the euclidean line 
through them. If L does not intersect the boundary circle C of 
the hemisphere II, then we take for E the vertical plane above 
L. 

Now suppose that L meets C, possibly tangentially. Note 
that P and Q are not separated by C as A and II are disjoint. 
Thus either both P and Q lie inside C, or both P and Q lie out- 
side C. Note that the second case has two subcases depending 
on whether P and Q lie in the same component of L — C. 
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Without loss of generality, we can suppose that \QC\ < \PC\, 
where the absolute value sign denotes euclidean distance. Let 
V denote the perpendicular bisector of PQ in the plane z = 
0, let O denote the centre of the circle C, let O' denote the 
intersection of V and the line OQ, let C denote the circle 
centred at O', through the point Q, and let £ be the hemisphere 
whose boundary is C . We leave it to the reader to check that, 
in all cases, E and IT are disjoint, as required. 
(7) Finding the distance between an algebraic geodesic segment e 
and a disjoint algebraic hemisphere II in hyperbolic 3-space can 
be algorithmically reduced to the preceding crises dbS follows. Let 
A denote the geodesic which contains e. If A and II are disjoint 
and do not have a common boundary point, we can apply 6) to 
find a lower bound for the distance between them. This is also 
a lower bound for the distance between e and II. Otherwise, 
the distance between e and II equals the distance between de 
and II, which reduces the problem to 4). 

Lemma 2.1. There is an algorithm to find a positive number R, such 
that star 11 (P 2 ) contains the metric ball B(P2,nR) . 

Proof. We will find R such that for each n > 1, the /^-neighborhood 
of dstar n (P 2 ), does not meet star n ~ 1 (P 2 ). Thus any path in H 3 which 
starts on dstar n (P 2 ) and ends on dstar n ~ 1 (P 2 ) must have length at 
least R. By induction it follows that any path which starts on dstar n {P 2 ) 
and ends on dP 2 must have length at least nR. It follows immediately 
that star n (P 2 ) contains the metric ball B(P2,nR), as required. 

We first give a description of the exact calculation before considering 
the error term. We use the above seven points to find positive lower 
bounds for various distances. The distance between disjoint vertices 
of P 2 can be bounded below using 1). The distance between a vertex 
v of P 2 and a disjoint edge e of P 2 can be estimated using 2), as v 
cannot lie in the geodesic which contains e. The distance between a 
vertex v of P 2 and a disjoint face F of P 2 can be estimated using 4), as 
v cannot lie in the plane which contains F . The distance between an 
edge e of P 2 and a disjoint face F of P 2 can be estimated using 7), as e 
must be disjoint from the plane which contains F. Finally the distance 
between disjoint faces E and F of P 2 can be estimated as follows. Let 
He and Up denote the planes which contain E and F respectively. If 
these planes are disjoint, and disjoint at infinity, we can estimate the 
distance between them using 5), and this will be a lower bound for 
the distance between E and F . Otherwise, the distance between E 
and F is bounded below by the distance between E and ITf. This last 
distance equals the distance between lip and some edge of E, and so 
equals one of the numbers already estimated. 

Now let R denote half the minimum of all these numbers. 
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Let W denote a vertex, edge or face of Pi. Then the definition of R 
implies that the ^-neighborhood of W meets only those vertices, edges 
or faces of P2 which meet W. 

Hence if Q is a translate of P2 in star n (P 2 ) which meets dstar n (P2), 
then the ^-neighborhood of dstar n (P2) does not meet any vertex, 
edge or face of Q except those which meet dstar n (P2). In particu- 
lar, it follows that the /^-neighborhood of dstar n (P2), does not meet 
star n ~ l (P2), as required. Note that dstar n (P2) and star n ~ 1 (P2) are 
disjoint. 

In the actual algorithm, when dealing with approximations, all the 
calculations above incur increasing error, but the fact that the number 
of operations is finite means that there is a constant C > such that 
the error in the estimate of each of these numbers is at most Ce, so 
that R must be replaced by R — Ce. It is of course possible that 
R < Ce in which case the algorithm must restart, replacing e by e/2, 
recalculating the coordinates in Manning's algorithm to this increased 
degree of accuracy, and recalculating R. Continue to do so until R > 
Ce, and then replace R by R — Ce, once this number is positive. □ 



3. Algorithm to find the geometry 

In this section we consider closed orientable irreducible geometric 3 
-manifolds, given by finite triangulations. There are eight geometries, 
and it is well known that a closed 3-manifold can have a geometric 
structure modelled on at most one of these geometries. This can be 
proved by exhibiting properties of the fundamental groups which dis- 
tinguish the geometries. For example only closed manifolds modelled 
on S 3 can have finite fundamental group. The point of what we do in 
this section is that we can decide algorithmically on which geometry a 
given geometric manifold is modelled. Such an algorithm is described 
by Manning in [Tlj 1.4], but here we provide a more algebraic treat- 
ment. 

As we are considering orientable irreducible 3-manifolds, the geom- 
etry S 2 x R cannot occur. For the only closed orientable manifolds 
modelled on this geometry are S 2 x S 1 and IRP 3 #RP 3 , neither of which 
is irreducible. We start by listing the remaining seven geometries to- 
gether with some selected properties of the closed manifolds modelled 
on these geometries. 
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Geometry 


Selected properties of any closed orientable 3-manifold 
M modelled on given geometry 


S 3 


7r%(M) is finite. 


E 3 


ni(M) is virtually Z 3 , and M is Haken. Any two-sided 
incompressible surface in M must be a torus. 


e 3 


7Ti(M) has no subgroup isomorphic to Z 2 . 




A// i 1 ^ P{ Spifprt finrp m~t£ipp witn hA/TiPrnnlip n^mp nrnifnlH 

M is Haken, and contains an embedded incompressible 
hyperbolic surface. 


Nil 


M is a Seifert fibre space with Euclidean base orbifold. 


Solv 


M is Haken. Any two-sided incompressible surface in M 
must be a torus. 


SL 2 R 


M is a Seifert fibre space with hyperbolic base orbifold. 
A two-sided incompressible surface in M must be a torus. 



A crucial fact for us is that if M is modelled on one of the above 
seven geometries, then 7Ti(M) is bi-automatic, except in the cases when 
the geometry is Nil or Solv (|U chapter 12] proves automaticity). In 
order to apply the theory of bi-automatic structures, we first need to 
be able to decide whether M is modelled on Nil or Solv. 

The key topological algorithm we will need is that of Jaco and Oertel 
[5] which decides whether a given triangulated 3-manifold M is Haken. 
Their paper also shows how to decide whether M has an incompressible 
surface which is a torus, and how to find such a torus. We will use these 
algorithms several times in what follows. 

As usual T denotes the 2-torus S l x S 1 . We also need to be able to 
decide whether M is homeomorphic to T x I. If M is irreducible, this 
is a special case of algorithm 9.7 of [9]. Note that the algorithms of 
Jaco and Oertel also find essential discs, and in this case cutting along 
a properly embedded disc gives a 3-ball, which can be recognised (by 
Rubinstein [T7] and Thompson [2U] , or by Perelman's solution of the 
Poincare conjecture). 

Lemma 3.1 (c/. Scott, Thm 5.5). If M is orientable and is obtained 
from T x I by gluing T x {0} to T x {1} by some homeomorphism h, 
then M is geometric and is modelled on one ofE, 3 , Nil or Solv. 

Proof. The action of h on Hi(T) = 1? is given by an integer 2x2 
matrix A. We consider the trace, tr(A), of A. If |ir(A)| < 2, or if 
A = ±J, then A, and hence h, must be periodic, so that M is modelled 
on E 3 . If |£r(A)| > 2, then A has distinct real eigenvalues, so that M 
is modelled on Solv. If |ir(>l)| = 2, then A has a repeated eigenvalue 
equal to ±1. So long as A ^ ±1, this implies that M is modelled on 
Nil. □ 
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Remark 3.2. Suppose that we have found an incompressible torus T 
in an orientable 3 -manifold M , using the algorithms of normal surface 
theory, and that we have checked that cutting M along T yields a man- 
ifold homeomorphic toT x I. In this situation, one can algorithmically 
calculate the action of h on -ffi(T) = Z 2 , and so can decide on which 
geometry M is modelled. 



We now discuss the geometries Nil and Solv in more detail, and 
describe an algorithm to decide which geometry occurs. 

If M is modelled on Solv, then M is a bundle over a 1 -dimensional 
orbifold with fibre the torus. Thus either M is a bundle over S 1 with 
fibre the torus, or M is double covered by such a manifold. 

If M is modelled on Nil, there are several cases. If M is Haken, 
then it is a Seifert fibre space whose base orbifold is a torus, Klein 
bottle, S 2 (2, 2, 2, 2) or P 2 (2, 2). If this orbifold is not a torus, there is a 
regular cover of M of degree 2 or 4 whose base orbifold is a torus. If M 
is not Haken, it is a Seifert fibre space whose base orbifold is S 2 (p, q, r), 
where (p,q,r) is one of (3,3,3), (2,2,4) or (2,3,6). Now the orbifold 
fundamental group of S 2 (p, q, r) is the triangle group A(p, q, r), and in 
these cases, A(p, q, r) has a homomorphism to Z 3 , Z 4 or Z 6 with kernel 
isomorphic to Z 2 . Thus there is a homomorphism of 7Ti(M) to {1}, Z 3 , 
Z4 or whose kernel determines a finite cover of M which is a circle 
bundle over the torus. If M is modelled on Nil and is a circle bundle 
over the torus, then any two-sided incompressible surface in M must 
be a vertical torus, and if we cut M along such a torus, the result will 
be homeomorphic to T x I. 

We can apply the preceding paragraph as follows. For any triangu- 
lated closed orientable irreducible 3-manifold M, we can check whether 
there is a homomorphism of 7r 1 (M) to {1}, Z 2 , Z 3 , Z 4 , Z 2 x Z 2 , or Z 6 
whose kernel determines a finite cover of M with infinite first homology 
group. If this does not occur, then M cannot be modelled on Nil or 
Solv. If this does occur, we can check whether the covering contains 
an incompressible torus, and if it does, we can check whether cutting 
along this torus yields T x I. If this occurs, the remark above tells us 
how to determine the geometry on this finite cover and hence on M. 
If it does not occur, then M cannot be modelled on Nil or Solv. 

Thus we can decide whether or not M is modelled on Nil or Solv, 
and if it is so modelled can decide which. 

This reduces us to considering the five remaining geometries. 

Given a triangulated closed orientable irreducible 3-manifold M, we 
can write down a presentation for %i(M). If we know that iri(M) is 
bi-automatic, we can algorithmically find a bi-automatic structure. We 
can then algorithmically check whether 7Ti(M) is finite. Thus we can 
decide whether M is modelled on S 3 . 
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This reduces us to the four remaining geometries, which are E 3 , H 3 , 
ffxl and SL^. 

If M is modelled on E 3 , there are several cases. In all cases, M 
is Haken. It is a Seifert fibre space whose base orbifold is a torus, 
Klein bottle, S 2 (2, 2, 2, 2), P 2 (2,2), or S 2 (p,q,r), where (p,q,r) is one 
of (3,3,3), (2,2,4) or (2,3,6). Thus as for Nil geometry, there is a 
homomorphism of tti(M) to {1}, Z2, Z3, Z4, Z2 x Z2, or Z@ whose 
kernel determines a finite cover of M whose base orbifold is the torus. 
But now this finite cover must be a 3-torus, and so have free abelian 
fundamental group. 

Thus to decide whether or not M is modelled on E 3 , we simply check 
whether there is a homomorphism of tt\{M) to Z 2 , Z 3 , Z4, Z 2 x Z 2 , or 
Z 6 whose kernel is free abelian of rank 3. For M is modelled on E 3 , if 
and only if there is such a homomorphism. 

This reduces us to the three remaining geometries, which are H 3 , 
l 2 xl and SL^. 

Using the bi-automatic structure on tti(M) one can check the an- 
swers to the following questions. Does m(M) have nontrivial centre? 
Does 7Ti(M) have a subgroup of index 2 with nontrivial centre? If the 
answer to both questions is negative, then M must be hyperbolic. If we 
find a positive answer, then M must be modelled on one of H 2 x R and 

SL 2 M.. To distinguish these cases, we use the facts that if M is modelled 

on SL2R, then any incompressible surface in M must be a (vertical) 
torus, whereas if M is modelled on H 2 x M, there must be horizontal 
incompressible surfaces in M none of which can be a torus. Thus we 
can apply the algorithm of Jaco and Oertel [5] to decide whether M 
contains an incompressible surface which is not a torus. 
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